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Abstract — Extremely high temperature achieved in the 
contact surface between the disc and the pad during braking 
has a large impact on the brake performance. Unintended 
consequence of extremely high temperature can cause complete 
brake malfunction, so called fading. There are a number of 
mathematical models describing the main tribological 
properties of various materials, including friction materials 
used as brake linings or pads. The application of these models 
for prediction or estimation of friction properties under 
determined working conditions is, however generally very 
limited. It is demonstrated in the paper by an example which 
has rather high practical importance, in order to shorten the 
experimental test. 

Index Terms — braking, temperature of friction surface, 
coefficient of friction, mathematical modeling. 


I. INTRODUCTION 

Numerous research showsthat the temperature achieved 
during braking on the contact surface between the disc and 
the pad has a large impact on the coefficient of friction [1-5]. 
Extremely high temperature, finally, can cause complete 
brake malfunction, so called fading [6]. A large number of 
authors used different mathematical methods in trying to 
describe a coefficient of friction in the contact zone of the 
friction pair, but a great majority of them are often very 
complex and contain a number of limitations [7-10]. The 
functional characteristics of friction materials for motor 
vehicle brake linings determine the brake performance as 
well as the behavior of motor vehicle during braking. It is 
also well known, that the best way to determine operational 
characteristics of friction materials is purely empirical that is 
to say experimentally. However, one knows that this method 
in the specific research area is lengthy and expensive. 

Therefore, theoretical and mathematical determination of the 
fundamental tribological characteristic of friction materials 
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in function tribological parameters, is of prime importance, 
due to fact that time and auxiliary expenses are radically 
reduced. In hitherto published works [11-14], sufficient 
space is given to the application of the mathematical 
modeling, predicting the operational characteristics of 
friction materials used for motor vehicle brake linings. This 
work has no intention to show explicitly the method of the 
actual mathematical modeling, but to direct itself on 
explaining how to develop the specific tribological model 
based on short term experimental results obtained from 
laboratory test of motor vehicle friction pads, on inertia 
dynamometer, as well as to explain the potential exploitation 
of the resulting tribological model in predicting the nature of 
friction under extremely high operating temperatures. 

II. Experimental research 

Experimental research was undertaken using frictional 
brake pad lining of motor vehicle, at Single-end full-scale 
inertial dynamometer PSK-20, developed at the Faculty for 
Mechanical Engineering, University of Belgrade. The 
outlook of the inertia dynamometer is shown in Fig.l. This 
particular inertia dynamometer is universal, i.e. it could be 
applied in testing the whole range of different brakes built in 
motor vehicles-cars, vans pick-ups etc. Single-end full-scale 
inertial dynamometer PSK-20 consists of three basic groups: 

1. the propulsory group, 

2. combination of different flywheels and 

3. the part in which the brake under the test is 
mounted. 

As shown in Figure 2 single-end full-scale inertial 
dynamometer PSK-20 is arranged ‘in-line'. A system of six 
flywheels (3) is mounted onto a common central shaft (4) to 
simulate the mass of the test vehicle. Variable inertia is 
obtained by engaging the appropriate flywheels, 10-200 
kgm . Precise inertia simulation is obtained by ‘trimming’ 
using the DC electric motor (1), with continuous regulation 
of angular speed which drives flywheels via a periflex clutch 
(2). The brake disc (8) is with flange (5) firmly attached to 
the shaft (4), while the brake calliper (9) is mounted for fixed 
flange (6) which is firmly connected to the foundation (7). 
The full-scale inertial dynamometer has been used for strictly 
controlling brake’s operating conditions namely application 
pressure, initial speed, and temperature at the contact 
surface. Whole system is equipped with a PC-based 
automatic control and data. 
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Figure 1 . Single-end full-scale inertial dynamometer 


Control unit manages angular speed and control line 
pressure, while data acquisition system measures: 
application pressure, speed, friction surface temperature, and 
braking torque. Brake torque is measured by strain gages, 
control line pressure by pressure transducers, speed by an 
optical encoder and temperature using thermocouple. 



Figure 2. Principal scheme of the measuring system for brake 

testing 

In this way repeatability of brake’s operation conditions 
for different types of friction material have been provided. 
The experimental research programmed consisted of two 
parts. The first part of the research comprehended the 
evaluation of the relationship between the friction coefficient 
(p) and activating pressure (p), speed (v) and temperature (0), 
so that the experimental results for the friction coefficient 
could be obtained, and used in the development of the 
specific friction model. In experiment were used 3 different 
samples brake pads (3 different producers) of same brake. 
Within the frame of experiments, 27 brake applications were 
performed, using different values of the influential 
parameters namely activating pressure, speed and 
temperature. The magnitudes of those parameters 
correspondent to the real braking conditions, experienced by 
the friction pad lining on a motor vehicle. The results are 
shown in table 1. As it can be seen in Table 1, measurements 
were carried out during 27 consecutive (or repeated) full- stop 
brake applications with an initial brake speed of 30, 60 and 
90 km/h, control line pressure of 20, 40 and 100 bar and 
initial brake temperature of 50, 100 and 150 °C. It is done in 
order to provide relatively average braking regimes under 
different brake’s friction surface temperature regimes. In 
each brake application, brake disc was first accelerated until 
reaching the predetermined initial brake speed and then 


braked to the full stop. After completing the braking, brake 
disc remains stand still. Speeding up of the disc for the next 
brake application starts few instances before the brake 
reaches predetermined initial temperature. 

Table 1. Results on friction tests 


Line 

pressure 

Initial 

brake 

speed 

Mean 

interface 

temperature 

Coefficient of friction 
EXPERIMENTAL (-) 

bar 

km/h 

°c 

Pad 

Pad 

Pad 




1 . 

2. 

3. 

20 

30 

50 

0,413 

0,425 

0,395 

40 

60 

50 

0,405 

0,431 

0,415 

100 

90 

50 

0,452 

0,441 

0,424 

20 

30 

100 

0,405 

0,420 

0,392 

40 

60 

100 

0,425 

0,435 

0,418 

100 

90 

100 

0,453 

0,438 

0,422 

20 

30 

150 

0,375 

0,421 

0,375 

40 

60 

150 

0,425 

0,436 

0,407 

100 

90 

150 

0,435 

0,445 

0,421 


In the second part of the experimental research, 63 brake 
applications were performed with the same activating 
pressures and speed, in the first part, while the operating 
temperatures were considerably above those normally 
developed in the process of braking on a vehicle under 
typical conditions (from 200 C to 500 °C). The results from 
this section of the research are shown in table 2. 


Table 2. Results on friction tests 


Line 

pressure 

Initial 

brake 

speed 

Mean 

interface 

temperature 

Coefficient of friction 
EXPERIMENTAL (-) 

bar 

km/h 

°C 

Pad 

Pad 

Pad 3. 




1 . 

2. 


20 

30 

200 

0,392 

0,415 

0,388 

40 

60 

200 

0,425 

0,419 

0,408 

100 

90 

200 

0,442 

0,444 

0,435 

20 

30 

250 

0,391 

0,405 

0,392 

40 

60 

250 

0,406 

0,420 

0,395 

100 

90 

250 

0,452 

0,429 

0,414 

20 

30 

300 

0,399 

0,402 

0,365 

40 

60 

300 

0,406 

0,422 

0,385 

100 

90 

300 

0,428 

0,435 

0,421 

20 

30 

350 

0,389 

0,386 

0,385 

40 

60 

350 

0,411 

0,415 

0,406 

100 

90 

350 

0,435 

0,445 

0,411 

20 

30 

400 

0,392 

0,403 

0,385 

40 

60 

400 

0,415 

0,418 

0,399 

100 

90 

400 

0,425 

0,425 

0,415 

20 

30 

450 

0,371 

0,376 

0,376 

40 

60 

450 

0,405 

0,413 

0,395 

100 

90 

450 

0,427 

0,441 

0,425 

20 

30 

500 

0,369 

0,404 

0,375 

40 

60 

500 

0,395 

0,411 

0,401 

100 

90 

500 

0,422 

0,438 

0,416 


III. MATHEMATICAL modelling ol lriction 

The aim of the first part of the mathematical modeling of 
friction was to find out the relationship between the friction 
coefficient (p) and activating pressure (p), speed (v) and 
temperature (0). To achieve this purpose the theory of 
experimental tests was used. The modern theory of 
experimental tests consists of the statistical multifactor 
analysis plans. These multifactor plans can be used for: 
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• the mathematical modeling of phenomena, 
processes and systems in space and time, 

• the study of the nature and mechanisms of internal 
processes, and 

• the optimization and optimal control processes in 
technical systems. 

Multifactor experimental design is characterized by two 
main features: 

• minimum number of experimental points, resulting 
in much lower costs and shorter duration of 
experimental investigations of processes and 
systems, and 

• maximum of information on the effects of a 
mathematical model of the process. 

Mathematical model of friction is realized using a 
statistical method known as regression analysis. Regression 
analysis deals with setting up a stochastic model of the object 
of research, which explains the condition and behavior of the 
given object, in a sufficiently reliable manner, within the 
covered experimental space. Regression analysis shows a 
form of connection between two variables using the 
regression line, which is determined by this condition: the 
sum of squares of the vertical distances of the points to the 
direction is the smallest, the so called method of the least 
squares. The least square method is an arithmetic procedure 
which allows adjusting the pre selected shape of functions to 
empirical distribution of the data and defining its parameters 
minimizing the sum of squared deviations of the functions of 
the empirical data sets. This method, however, does not 
reveal the shape of the function which is best suited to 
empirical data, but allows a choice between the functions of a 
given shape which suits best. Results obtained in any 
experimental research allow to determine interrelation 
between the test results in mathematical form. In planning of 
the experiment, there are two cases: 

• mathematical model of tests phenomenon or system 
is known and based on previous knowledge of the 
limited experimental area, or 

• mathematical model is unknown 

The actual analytical form of response function is actually 
unknown, and the mathematical form is a more or less 
accurate model of it. During the research, a numerous 
mathematical models of friction are formed. The coefficient 
of friction during braking of motor vehicles, is variable and 
depends on many factors, which is well known. Basically, 
influential factors on the functional characteristics of friction 
materials for brakes of motor vehicles can be classified into 
two groups. The first group consists of factors that depend on 
the operating conditions or methods of use of the brake 
system (pressure, temperature and speed), and the second 
group consists of construction impacts, friction material and 
metal elements. With this in mind it is assumed that the 
appropriate mathematical model of the friction coefficient 
applied can be a stepped function of the third order in these 
forms: 

Ml =k r p A ^v B ^6 c ^ (1) 

ju 2 — k 2 * p A2 • v Bl • 0 Cl ( 2 ) 

jU 3 =k 3 -p A i -v* 3 (3) 

Where: 

p - activating pressure, 

v - speed, 

0 - temperature, 


k b k 2 , i k 3 - constants that depend on the structure of the 
applied friction material and applied metal element and 
Ai s A 2 , A 3 , Bi 5 B 2 B 3 , Q C 2 i C 3 coefficients to be 
determined experimentally. 

The assumed function of the coefficient of friction must be 
expressed logarithmically 

M,2,3 = A, 2 , 3 A, 2 , 3 l* 1 P + #1,2,3 ^ V + ^1,2,3 ^ 9 ^ 

and apply a method of the least squares to obtain constants 
k 123 , and coefficients A 123 , B 123 and Ci >2>3 . Using the 
experimental result for the friction coefficient obtained in the 
first part of the research, the specific method of mathematical 
modeling was conducted, with the following resulting 
expression: 

(5) 

( 6 ) 


Pad 1.: M =0,372 -/A 068 . v ^ 017 . < 9^0202 


Pad 2.: =0,384*/? 

Pad 3.: p 3 =0,365*/? 


0,07 - 0,02 

* V 

0,063 , ,- 0,016 


V 


0 


0 


- 0,195 


- 0,0178 


(7) 


IV. ANALYSIS OF OBTAINED RESULTS 

The fidelity of such friction model is illustrated in fig. 3, 
where one can clearly see that the differences between the 
friction coefficients obtained theoretically and 
experimentally differ quite negligibly. Figure 3 shows the 
results comparison of experimental and theoretical values of 
the coefficient of friction on the first sample of pads. In all 
27 brake applications differences were below 3%. It is also 
important to note that the friction model shown by the 
expression (1, 2, 3) characterizes the behavior on the friction 

coefficient in function of the influential parameters (p, v, 9 ). 





Figure 3. Comparison of experimental and theoretical values 

of the coefficient of friction 
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After verification of the mathematical model of friction by 
comparing theoretical and experimental results carried out is 
estimation the value of the coefficient of friction under 
extremely high temperature. By comparing the results 
obtained in the second part of the experimental research, 
shown in table 2, and the results using expression (1,2, 3), it 
can be seen that the actual error is less than 3%, which is 
satisfactory. This can be clearly seen on fig. 4. 




Q 



Figure 4 shows the comparative values obtained by 
experimental tests and the results obtained by mathematical 
model for the first brake pad. At this stage, it is quite easy to 
realize that the values of the friction coefficient, resulting 
from equations 5, 6 and 7 and experimentally ( at increased 
temperatures), are, by its character, in function of the 

influential parameters (p, v, 6 ) equivalent, which suggests 
that the friction model resulted from short term laboratory 
testing on inertia dynamometer, is virtually applicable for 
predicting the friction of the applied friction material, even at 
extremely high temperatures, on the same vehicle, under 
quite unfavorable operating conditions. 

V. CONCLUSION 

In the first part of this paper results of the experimental 
tests were used to develop a mathematical model of friction. 
Using the mathematical model of friction, value of the 
coefficient of friction under conditions of extremely high 
temperatures it was estimated and that was confirmed by 
experimental results in the same conditions. On the basis of 
conducted experimental research, obtained results and their 
analysis, the following conclusions can be drawn: 


• The results from 27 brake applications enable the 
developing of the friction model for the given brake pad 
lining. 

• Reliability of the results presented model is confirmed 
by experimental results which show the deviation is in 
frame 3%. 

• The resulting friction model can serve in predicting the 
friction and in that way the efficiency of the brake pad 
lining material, even at increased temperatures. In this 
way estimation of the actual behavior of the automotive 
brake is ensured in exploatation, even under extremely 
difficult operating conditions. 
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